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ABSTRACT 



A pilot engaged in a strike mission against targets of opportunity 
is continually confronted with the decision as to how many and which of 
his available weapons should be expended against a variety of targets 
which are encountered*, Associated with this decision problem is the 
problem of determining an optimal load for such a mission subject to 
the constraints of available payload and number of weapons stations on 
the aircraft® Certain assumptions are made concerning the distribution 
of targets within the target area which lead to a dynamic programming 
formulation of the decision problem® This yields a system of ordinary 
differential equations which are solvable recursively® In addition to 
the dynamic programming model , a sub -optimum determination of a "best 
load" is discussed® Although not as complete or precise as the dynamic 
programming method*, this formulation is more readily adaptable to squa- 
dron level decision making® 
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1. INTRODUCTION 



The problem of selecting weapons for an aircraft strike mission 
and using those weapons effectively has been documented quite exten- 
sively. However, present publications are devoted almost exclusively 
to the situation where the target is assumed known prior to the 
mission and the aircraft can be loaded for that target, or where 
several different aircraft can be loaded differently to cover an 
unknown spectrum of targets."^ These approaches are inadequate for 
determining weapons requirements against targets of opportunity. A 
pilot engaged in such a mission is confronted with the decision as 
to which and how many of his available weapons to expend when any 
target is encountered. When various types of targets might be en- 
countered, each with a different strategic or tactical value, he 
wishes to maximize the value of the targets he can expect to destroy. 
Associated with this is the problem involving the choice of weapons 
that should be initially loaded onto the aircraft for such a mission. 
Here the problem is constrained by maximum payload, weapons, and wea- 
pons stations availability. It is desired to obtain, within these 
constraints, (1) the optimum feasible weapons mix to be loaded for 
use against an uncertain spectrum of targets, and (2) a doctrine for 
the expenditure of these weapons when targets are encountered. 

The model developed herein attacks these problems in reverse 
order. First, an optimum firing doctrine is generated for all feas- 
ible loads. Then that load which yields the maximum expected return, 

^NAVWEPS Report 8664, "Preferred Aircraft Loadings Based on the 
Predicted Combat Effectiveness of Conventional Air- to-Surface Weapons," 
NOTS TP 3710, 1965, (SECRET). 
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when expended in accordance with the optimum doctrine, becomes the 
optimal load* 



ASSUMPTIONS IN THE MODELS 

Certain assumptions, necessary to the model development are made 
concerning target distributions within the target area and tactics of 
the strike aircraft. 

It is assumed that there are a finite number of different target 
types, each randomly and independently distributed within the target 
area with known density. For example, from photo-reconnaissance or 
other sources it might be known that the density of trucks within the 
target area is x per square mile. The number of targets of each type 
is assumed to be large enough so that the probability of future en- 
counters is unchanged by consideration of the number already met. 
Further, these targets are assumed to be of a unitary nature, charac- 
terized by the existence of single shot kill probabilities, as dis- 
tinguished from area targets characterized by probabilities of des- 
troying a certain percentage of the area or value of the target. 

The assumption is made that a ratio scale of target values can 
be specified. For example, with the least valuable target as a base 
with value V, the achievement of a kill on a target type i is worth 
7V, target type j is worth 8V, etc. 

For defensive reasons the strike aircraft is restricted to a 
single attack on any target which might be encountered. The pilot 
may not drop a bomb on the target, assess damage, and then make 
another attack if a kill has not been achieved. He may however, 
expend any number of weapons on this single pass provided they are all 
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the same type. That is, while he may expend three 500 lb. bombs, he may 
not drop one 500 lb. bomb and one 1,000 lb. bomb. This restriction 
is not necessary to the development of the models, but rather recog- 
nizes that differences in aim points between different weapons make 
this mixed drop undesirable. 

METHOD OF SOLUTION 

Using data on weapons weights, constraints of aircraft payload 
capacity, and the number of available weapons stations on the plane, 
all feasible loadings can be enumerated. Initial loads are then con- 
sidered which are not subsets of some other feasible load. That is, 
if five weapons of type one and three weapons of type two constitute 
a feasible load, then four weapons of type one and three of type two 
is also feasible but a subset of the former. The kill potential of 
some subset of a feasible load is clearly less than that of the 
original load and is therefore not considered in the determination 
of an optimal initial loading. 

A functional equation for the expected value of targets killed 
is then set up as a dynamic programming process yielding a value for 
kill potential for each of the feasible loads being considered. The 
maximum kill potential from amongst these loadings determines the 
optimal load for the strike aircraft. Additionally in the solution 
of this dynamic programming process, the optimal firing strategy for 
the pilot is obtained for any initial load. The optimized sequential 
decision process yields a policy which determines as a function of 
the operating time and weapons remaining, the best weapons to be ex- 
pended in an encounter with any of the various target types. 
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Dynamic programming provides a complete and very precise means 
of determining the optimum load for a "targets of opportunity 11 type 
mission. Except in the very simplest cases however, the problem is 
quite complex and will require a computer for solution. In the majority 
of operational situations such equipment will not be available. Thus 
an approximate method, suitable for hand calculation, is needed for 
determining the initial load. 

In order to achieve simplicity, this "hand model" must sacrifice 
detail and thus completeness. In this respect it is a sub-optimal 
method. The simplified method is best suited for comparing a few, 
rather than all, feasible loads for use against a selected target 
list. By choosing the best weapon for each target type and comparing 
the expected returns for various mixes of these weapons, a "best load" 
can be generated. 
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2. GENERAL DISCUSSION AND DEVELOPMENT OF THE MODELS 



The random distribution of targets leads to an exponential inter- 
encounter time distribution for the strike aircraft searching the 
target area.^ (All targets which are within the visibility range of 
the aircrafts track are assumed sighted. Alternatively the sightings 
could be reduced by some uc 0 ^auation factor which takes into account 
the difficulty associated with perceiving a given target type. This 
however would not alter the analysis which follows.) The mean rate 
of sightings (the reciprocal of the mean of the exponential distri- 
bution of inter-encounter times) is denoted by r. 

DYNAMIC PROGRAMMING MODEL 

The development of the dynamic programming formulation of the 

problem is best seen by first considering the simpler cases. To begin 

with, a model has previously been developed for the allocation of wea- 

2 

pons over a limited time interval. This can be utilized for gener- 

ating the optimal firing strategy in the simplest case of one target 
type-one weapon type* The optimal loading problem is not pertinent; 
for a single weapon type the best load is clearly the maximum number 
that can be carried on the aircraft. First define: 

f (t) = expected number of targets killed in the time 
n interval (0,t) if there are n weapons at time 

t = 0 and these are utilized optimally. 



1 

Appendix I 

2 

Bram, J. , "Allocation of Weapons to Targets with Exponential 
Arrival Times in a Limited Time Interval, 1 OEG Interim Research 
Memorandum 32, 1963. 
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An initial loading of N weapons and a time T allowable for search 

is assumed, and f (t) will be calculated for 0 ^ n ^ N and 0 ^ t ^ T. 
n 

Define: 



p « probability that j weapons expended against a target 
j will achieve a kill. 



Clearly f (t) is dependent on the p 1 s and on r, the encounter 

n j 



rate*. 



By Bellmans "principle of optimality,"^ for 0 ^ t^ t, the 
optimal policy in the interval (0, t) , given n weapons at time t = 0, 
must also be optimal in the sub interval (t^, t) given n - n^ weapons 
at time t = t^, where n^ weapons were expended in the first subin- 
terval (0, t^) • 



Let t^ = At « T. By definition of the exponential arrival pro- 
cess, the probability of an encounter in a small interval of time At 
is proportional to At with constant of proportionality r, the en- 
counter rate. Thus the probability of one encounter in the time 
interval (0, At) is rAt and no encounters occur with probability 
1 - rAt. 

Breaking the time interval (0, t) into the two sub- intervals 
(0, At) and (At, t) yields: 



(1 - rAt) = the probability that no encounters will occur 
in the first sub-interval. 

p. = the expected return in the first sub-interval 
J having expended j weapons. 



Bellman, R„, Dynamic Programming , (Princeton, New Jersey: 
Princeton University Press, 1957), p. 83. 
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f (t - At) « the expected return in the second sub-interval 



having expended no weapons in the first and following 
an optimal strategy in (At, t) . 



f (t - At) ° expected return in the second time interval having 
n ~j expended j of the n weapons in the first sub- 

interval and behaving optimally in (At, t) . 



The expected return over the whole interval (0, t) is then: 



To optimize over the period (0, t) , j must be chosen such that this 
expression is a maximum. Therefore 



The maximization here is over the range l^j^n rather than O^j^n 
since it clearly would never pay to let a target go by without firing 
any weapons, saving them all for future encounters which are not 
certain to occur. 

Rearranging (1), dividing by At and taking the limit as At 
approaches zero yields: 



From (2), f (t) can be calculated for all n and t recursively 
n 

as follows: By definition, f (t) = 0 for all t. (The expected 

0 

return is zero over any period of time if there are no weapons with 
which to operate.) Also f n (0)= 0 for all n. (The expected return 
is zero if there is no time in which to operate.) Solving (2) for 



(1 - rAt) f (t - At) + rAt (p + f (t - At)) 
n j n-j 



( 1 ) 



f (t) = (1 - rAt) f (t - At) + rAt max (p + f (t - At)) 
n n j n-j 



l£ j^n 



( 2 ) 
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n = 1 yields: 



f x (t) = Pl (1 - e‘ rt ) 

Note that this result was to be expected. Exponential inter-encounter 
times lead to a Poisson counting process for the number of encounters 
within a time interval. The term (1 - e c ) , is simply the prob- 
ability of at least one encounter in the time interval (0, t) , and p^ 
is the probability that, given an encounter, a kill will be achieved 
with the expenditure of the single weapon. 

Utilizing f^(t) equation (2) can now be solved for n = 2. 



f'(t) = - rf (t) + r max (p (2 - e rt ) , p ) 

2 2 1 2 

-ri- 1 P 1 

Since p^> p^(2 - e ) when either t < ~ In or if p^> 2p^, 



the solution to the differential equation is: 

p l 



f (t) ■ P 0 (1 - e" rt ) for t ^ “ In 
2 2 r 



or if p s 2p, 

2p r p 2 2 1 



f 2 (t) = 2 p l - e'" (r Pl t + Pj + P 2 - Pl In 2p ^ ) 



for t s — In 
r 



2p -p 
1 2 



and p 2 < 2p x 



The significance in this solution lies in the "cross-over" point in 
time. If the search time remaining is less than this value of t, the 
optimal return is achieved by dropping both bombs on the next target 
which is encountered. If greater, only one bomb is dropped on the 
target, the other being saved for possible future encounters. 
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Having found f^(t) , f^(t) can now be calculated. Tlie differential 

equation however becomes unduly complicated for n > 2. A numerical 

solution may be obtained for the problem though by recourse to equation 

(1), the difference equation, with which f (t) can be calculated for 

n 

any n to any desired degree of accuracy by taking At small enough. 

With the possibility of encounters with targets of different 
types, the value of a target which is met must be considered when the 
decision is made as to how many weapons should be expended. 

The methods used in developing the functional equation are 
analagous to those for the simpler model with the addition of con- 
siderations of target value and different kill probabilities. Define: 



m » the number of different target types. 

p =» the probability that j weapons expended on a target 
type i will achieve a kill. 

r, * average encounter rate with type i targets, 
v^ « value of target type i. 

f (t) * expected value of targets killed in the interval (0, t) 
given n bombs at time t = 0 and optimal utilization of 
these weapons. 



Then: 



(3) 



m 

f (t) = (1 - r At) 
n I— x 

i=l 



f (t - At) + . r At max (v^p +f ,(t-At)) 
n ' U x 1 ij n-j 

i=l O^j^n 



Note that in this case the maximization is over the range 0 £ j ^ n 
vice 1 ^ j ^ n. Here the expected return might be greater if no 
weapons are dropped on targets of little value, expecting at some 
later time in the mission to meet something of greater worth. 
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The differential equation is obtained by subtraction of f^(t) 
from both sides of the above equation, division by At and taking the 

limit as At approaches zero* 

m m 

V 1 V 

f'(t') = -) r f (tl + } r. max (v.p. . + f .(t)) 

V ' U i + u i _ ^ i ij n- j v 

i= i i=i ° 5 J^ n 



For n = 1 and m = 2 the solution to this equation can be shown tc be: 



f l (t) 



and 



< r lVl + VaV n -(r + r_) t . 

(1 - e 1 2 ) for p v > p v 

(r. + r_) 2 2 1 1 



O^t^ 



1 2 
1 



r i ' ! ' r 2 



ln (r l p l v l + r 2 p 2 v ?> 
” r 2 <P 2 V 2 ' P 1 V 1 > 



or if p v,> p v then for 
1 1 2 2 



1 , (r i P l v l + r 2 P 2 V 2 ) 

0£t£ in =-=— 



I - + r„ , y. 

1 2 r iVl • P 2V 



or if p v = p v then for all t. 
11 2 2 



£ 1 (c) = p v + (p v - p v > [ ^ P -- 2^2_2 1 r 

11 22 11 SVi-Yj* 



+r e 
2 



'V 



r p,v, + r p v 
1 111 2^2 2 
for p v > p v and t ^ — T — in 

1122 r l r 2 r !^ p l V l ’ P 2 V 2^ 



f n (t) = p v + (p v -P v ) , - — " 

1 22 11 22Lr 2 (p^ - PjVj) 



r r i p i v i + r 2 p 2 v 2 j 7^1 - r 2 t 



1 r l p l v l + 9 P ? V 2 

for p v > p-v. and t £ — - — In — *■ — 

2 2 *11 r +r r (p v • p v.) 

1 2 2 2 2 1 1 
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Thus for the simple case of two target types and a single weapon, 

the expected returns can be calculated analytically. In addition, the 

optimal strategy is found, here consisting merely of a cross-over time, 

the time before which it does not pay to fire the weapon at the target 

with the lowest p v . For n > 1 however, it can readily be seen that 
i i 

solutions to the differential equation are not so easily obtained. 
Again though, through recourse to the difference equation, (3) , numer- 
ical solutions can be obtained for any n to any desired degree of 
accuracy. 

Further generalization of the model entails a relaxation of the 
restriction to a single weapon type. The functional equation is 
similar to those for the simpler models with some minor changes in 
the definitions of the terms. Define: 



p = the probability that k weapons of type j expended on 
ijk a target of type i will achieve a kill. 

L * the number of different weapon types. 

f (t) = expected value of targets killed in the Interval 

n ^ n 2 n ;j** #n L (0, t) , given n. weapons of type i (for i"l,»**,L) 

at t = 0 and optimal utilization of these weapons. 



The difference equation can be written: 



f 

n n n . 
12 3 




m 

= (1 r At) 

L - l 

i=l 



f (t - At) 

n n n . . . n 
12 3 L 




r At 
i 



i=l 



max max ( V .P, + f 

i ljk n n 

l£j£L O^k^n 1 2 

j 



(t - At)) 

. n • • n 
j L 
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The differential equation can be written as before, but the 



solutions are obtainable only in the most trivial cases • 

A FORTRAN program has been written however from which numerical 

solutions to the difference equation can be calculated."^ With inputs 

of the aircraft load carrying capabilities, the mission time, and the 

parameters of the model (m, r., v , L, p , and At) the program 

i i ljk 

selects values for the n which are consistent with the limitations 

i 

of the aircraft and maximize f . Additionally, another FORTRAN pro- 
gram'*' lists the optimal strategy to be employed in the decisions 
relating to the expenditure of these weapons. 



A SUB-OPTIMAL DETERMINATION OF LOAD 

In the majority of practical cases hand calculation will be the 
only means available for determining loading requirements. In such 
cases a sub-optimal method can be used at the squadron level for 
obtaining the initial weapons load for a mission* Define: 



T = search time (i.e., mission time allotted to searching for 
and attacking targets). 



P . (T) 
ni 



probability that target type i will be met n times during 
the mission. 



n ,(T) =^P .(T) 

1 ' r 



n=l 



ni 



probability that target type i will be met at 
least once during the mission. 



v “ value of target type i. 
i 

ra = number of target types. 

b_. = number of sticks^ of s weapons of type j in the load. 

L = number of weapon types available. 

^*In this development the terms stick and salvo are used equiva- 
lently. The choice of a stick or a salvo, in the load determination, 
will depend upon which gives the higher kill probability. 
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probability that a stick of s type w weapons will kill 
target type i. 



k 

iws 



Note that in this development individual bombs are considered as 
a stick of one; i.e., two 250 lb, bombs might be considered as two 
sticks of one type j or one stick of two type j weapons. 

If, a priori, only the more highly valued targets and those 
individual weapons with the highest kill probabilities against these 
targets are considered the assumption can be made that a target 
sighted will be dropped on, as long as the appropriate weapon for 
that target is available. (Note that this assumption is used only 
for determining the load. In the actual performance of the mission 
there is the possibility that, if that weapon most effective against 
a target is not available, a less effective weapon will be dropped.) 
The problem then becomes one of maximizing, for a given mission time, 
the expected value of targets killed, E, for those loads being con- 
sidered. E may be determined as follows: 
m 

p D f k v 

E L iJ s ijs i 
i-1 



where f is the expected number of type i targets on which weapon 
ijs 

type j can be expended, and j indicates that weapon type for which 

kj u £ k for all w = 1, 2, -## , L. 
iwl 

If the number of type j weapons available were limitless, f , . 

co ijs 

would simply be equal to n = ^ np^CD ' However > since a given load 

n=o 

is limited to b sticks of type j weapon, no more than b. can be 
js js 

dropped regardless of how many type i targets are sighted. Thus 



W decision as to which targets and weapons will be considered 
will depend highly upon the judgement of the individual responsible for 
determining the aircraft load. 
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f ijs is S iv en by 



(4) 



b . “ 1 oo 

JS 

f,,- " / np .(T) +Y b p (T) 
ijs L ni ^ js ni 

n=o n=b 

js 

Equation (1) can be reduced to a more easily calculable form. 



Since b . -1 b. -1 

D js js 

v 

, n P . (T) 

ni n 

n=o n=i 



b. "I 



-l ■“>„!<« ° nd I P ni <T> '"I' L p ni (1 >' 



n=b 



js 



n=l 



equation (4) can be written as 
"js - 1 

£ iis - I np (T) + b. n 
L JS L-> ni JS 1 

n=l 



- b 



k js”* 



V 



JS 



P .(T) 
ni 



n 



=1 



b. -1 
js 

= b n - £ (b . -n) p (T) 

js i “ J s ni 
n-i 



but II. = 1 - p (T) and consequently b II = b - b. p (T) . Note 
1 Oi js i js js Oi 



that b, p can be expressed as (b. - 0) p 

js Oi Js Oi 



Thus equation (4) 



now reduces to , , 

b - i 

js 

(5) f = b - ' (b. -n) p (T) 

ijs J s J s ni 

n=o 



Equations (4) and (5) are valid for all b ^ 1. By definition 

js 

f g = 0 for bj s = 0, since no sticks of type j weapon can be dropped 
if none are carried. 

As noted before, the random distribution of targets and conse- 
quent exponential inter- encounter times lead to a Poisson counting 
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process in determining the probabilities of encountering given numbers 

of targets; i.e., for type i targets appearing at a mean rate r , 

i 

n 



p .<T) - e' r i T 



ni 



_ i 

n. 



This gives 

( 6 > 



b -1 

js 



f = b 
ijs J s 



-r T frjT) 
(bi -n) e i 
J s n. 



n=o 



Table 11, Appendix III gives values of f_ $ for various com- 
binations of b and r T. This table greatly facilitates calculation 
js i 

of a "best load" from the feasible loads being considered. 

To determine the "best load" for a given set of targets, find 
that individual weapon type j, for each target, which has the highest 
kill probability. The loads considered are then the different mixes 
of these weapons selected from a list of feasible loads for the air- 
craft. 

For a given load, b. will be known for each type weapon in the 

J s 

load and the encounter rate r will be estimated for each of the 

i 

critical targets. The planned search time T will also be known for 

the mission. Thus using known values of k, , and estimated values 

ijs 

of v^, E can be calculated. 

Consider first the case where weapons will be dropped indivi- 
dually (sticks of one only) and where each weapon type in the load 
will be used against only one type target. In other words, target 
types and weapon types are matched one to one. Since weapons are to 

be dropped singly, b = b . 1 for each weapon type j. To determine E 

js j i 

for such a load, first determine b and r.T for each weapon- target 

jl 1 

type combination. Then obtain f from Table 11, and multiply this 

ijs 
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value by k v 
ijs i 



i 



to 



. Let f k v 
ijs ijs i 



E • Add the resulting E 

ij 



get E for the load, 



i.e. 



m 

"I'u 

i=l 



(Note that the E..*s are summed 
iJ 



over i only, since the subscript j simply indicates that weapon which 
is most effective against target type i.) In the same manner, E may 
be determined for each load being considered; that load which yields 
the maximum E is then the "best load. 11 

For this sub-optimal model, the extension of the problem to 
include sticks of two or more weapons is accomplished quite easily, 
since different size salvos or sticks of a given type weapon con- 
stitute different loads. For example, as noted previously, two 250 
lb. bombs might be considered as a load consisting of two sticks of 
one bomb or as another load consisting of one stick of two bombs. 

With this in mind, E can be calculated for these loads in the same 
manner as previously outlined. In order to keep the problem hand 
workable, this approach assumes that each delivery of a given type 
weapon will consist of the same size stick. This assumption is made 
solely for determination of the aircraft loading prior to flying the 
mission. During the actual mission, however, it is quite likely 
that sticks of different sizes will be expended at different times. 

For example, a pilot might drop one type j bomb on the first truck 
he sees, a stick of three type j bombs on the second truck he sees, 
etc. Prior to the mission, however, this load of four type i weapons 
could have been evaluated as one stick of four, two sticks of two, or 
possibly as four sticks of one, since calculations for sticks of 
varying sizes quickly exceed the bounds of practical hand workability. 
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In addition to the case where targets and weapons are matched one 
to one, consideration must be given to the case where a given weapon 
type is optimal for more than one target type. If the assumption that 
a target sighted will be dropped on were adhered to, the expected 
number of type i targets killed would be proportional to the encounter 
rate r^. However, in this situation, a concession is made to the fact 
that, early in the mission, a target of lesser value might be passed 
in hopes of encountering a more highly valued target later. Thus it 
is assumed that the expected number of type i targets killed is pro- 
portional to both the encounter rate, r., and the target value, v . 

i i 

The expected return from a given load is now calculated exactly as 

before with the following minor modification: For all targets type i 

against which weapon type j is most effective, calculate the E ^ j 1 s 

exactly as in the simpler case (as if all type j weapons were to be 

expended on each type target). Then to approximate the particular 

expected value of type y targets killed by weapon type j, E , is mul- 

r v v v 

tiplied by the normalizing factor N = — * — where the summation 

y s r i v i 

in the denominator is taken over those type i targets for which weapon 
type j is optimal. These expected values are then added to give E as 
before. 

A logical extension to this hand approximation is the consider- 
ation of "secondary" weapons. A secondary weapon would be considered 
to be a weapon which has high kill probabilities against several 
targets, yet is optimal against none of them. There are certainly 
cases where such a weapon might give higher expected returns over the 
entire target spectrum than would a load consisting of only those 
weapons which are optimal for the individual targets. For example, 
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consider a case in which there are two targets and three weapons to 
be used against these targets. The single stick kill probabilities 
for the individual weapons are as follows: 



WEAPON 


TARGET 1 


TARGET 2 


1 


• 

o 

vo 


.01 


2 


rH 

o 

• 


.09 


3 


.08 


.08 



If the load were limited to those weapons which are optimal for 
individual targets, only weapons 1 and 2 would be considered. However, 
it is obvious that weapon 3 merits examination, and, in fact, it appears 
as if a pure load of weapon type 3 might possibly be the best load to 
use in this particular case. Calculation of E for such a secondary 
weapon is identical to that in the foregoing case in which one weapon 
is considered for use against more than one target. 

This sub-optimal development may, at first glance, appear to be 
long and cumbersome. It must be realized though that in most actual 
situations all but a very few feasible loads will be eliminated from 
consideration prior to any calculation, and these will be examined 
only as far as they pertain to a small number of critical target types 
(normally not more than three or four). Thus it can be seen that 
relatively few calculations will actually be required and these will 
be of a trivial nature. 

In order to clarify the use of the sub-optimal model, the follow- 
ing simple example is presented in detail. 

Consider an armed reconnaissance mission with a planned search 
time of one hour. (Thus T = 1 so that r^T = r^ for all targets). 
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Three target types are considered: 



v 

1 



= 2 



= 5 



v = 1 
2 



r =3 
2 



v = 4 
3 



r 

3 



= 2 



The aircraft has two weapon stations and two types of weapons 
are considered, one more effective against target type 1, the other 
more effective against types 2 and 3 (see Table 2). 

The possible loads are indicated in Table 1. 



“1 

i 

1 




WEAPON | 

1 2 


Load 


Single 


Stick of 2 Single Stick of 2 


A 


1 


1 


B 


2 




C 


1 


2 


D 


1 


1 | 


E 




1 



TABLE 1 

COMPOSITION OF POSSIBLE LOADS 



, j 1 s - 


i 1 


2 


[ 

3 


1 
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.02 


1 .03 


.01 


1 ! 

2 


.06 


! .08 


.03 


1 

2 , 


.01 


.05 


.02 


2 


.03 


.12 


.05 



TABLE 2 

KILL PROBABILITIES 



27 



Load A: b 



1 



11 



= 1, b 



21 



Target 1 - Weapon 1 

f = .99326 (from Table 11; b. = 1, r T = 5) 

111 J s i 

E - f k v - (.99326) (.02) (2) = .0397 

11 111 111 1 

Target 2 - Weapon 2 

f = .95021 (b = 1, r T = 3) 

221 js i 

E = f k v = (.95021) (.05) (1) = .0475 

22 221 221 2 

Target 3 - Weapon 2 

f = .86466 (b = 1, r T - 2 
321 js i 

E = f k v = (.86466) (.02) (4) = .0692 

32 321 321 3 

Since weapon type 2 is optimal for two target types, E and E 



N E = (r v /rv + r v ) E = (3/11) (.0475) = .0129 
2 22 2222 33 22 



N E = (8/11) (.0692) = .0608 
3 32 

Thus for load A, E = .0397 + .0129 + .0608 = .1134 
Load B: b =2 



must be normalized 



11 



f = 1.95283 (b. =2, r T = 5) 

111 js i 



E 



11 



= E (Load B) = (1.95283) (.02) (2) = .0724 



Load C: b =2 

21 



r 



221 



- 1.75106 (b. = 2, r T = 3) 

J5 i 



i 



28 



N 2 E 21 “ < 3/11 > (1-75106) (.05) (1) = .0240 

f - 1.45866 (b. =2, r.T = 2) 

321 J s L 

N E = (8/11) (1.45866) (.02) (4) = .0846 
3 32 

E(Load C) «* .0240 + .0846 = .1086 

Load D: b “1 (Note that this is almost the same as Load B, but 

12 

now the two weapons are considered as one stick of two vice two sticks 
of one) . 

f = .99326 

112 

E - E(Load D) = (.99326) (.06) (2) = .1192 

Load E: b =1 

22 

f = .95021 

222 

N E « (3/H) (.95021) (.12) (1) = .0312 

2 22 

f = .86466 

322 

N E = (8/11) (.86466) (.05) (4) =-.1252 

3 22 

E (Load E) = .0312 + .1252 = .1564 

For the loads considered, the maximum E is .1564. Therefore the 
"best load" is load E which consists of one stick of two type 2 weapons. 
Note that if sticks of one only (case 1 in the model development) had 
been considered, the "best load" would have been load A with one of 
each type weapon yielding E = .1134. 

The foregoing sub-optimal method is limited in that it considers 
search time as a whole rather than a time-stepped firing doctrine as 
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in the dynamic programming method • Thus those added values which 
might be generated by strategies which vary with mission time are 
lost in this formulation. Consequently the expected kill values 
will generally be much lower than those of the dynamic program formu- 
lation, and the "best load" generated will not always be the optimum 
load. However, the load obtained by this method should, if expended 
in accordance with the optimum firing doctrine of the dynamic program, 
be one of the more highly valued loads. 
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3. sample results 



OPTIMAL LOAD AND F IKING STRATEGY DETERMINATION 



For a mission time of 15 minutes, seven weapons from which to 



choose, and values of the other input parameters as listed in Table 
3, an optimal loading has been determined using the ,, BESTLD m program 



of Appendix II, with the following inputs: 



Weapon Type : 1 2 3 4 5 6 

Weapon Weight: 265 300 510 600 675 950 

Number of racks “ 5 
Max payload ■ 5,000 lbs. 

Number of target types ■ 4 
Mission time - 0.25 hours 



7 

1975 



i 



Delta T ■ .010 houi‘6 (This is the At of equation (3)). 



TARGET TYPE 
1 

2 

a 

4 



ENCOUNTER RATE 
30 . 0OC0G 
1 . 000.0 
. 1 0< n.O 
3. Of 0'»0 



VALUE 
1 .0C 
2 5.00 
6 0 . 00 
5.00 



KILL PROBABILITIES 

1 O J 
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TARGET 
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WEAPON 
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. 0 2 C 
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WEAPON 
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. 049 


.06 0 
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.0 80 


TARGET 


I 


W E A PON 
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.040 


. 0 70 


.09 0 


. 100 


.110 


TARGET 


I 


WEAPON 
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. 00 5 


.008 


.01 0 


.011 


.012 


TARGET 


I 


W E A PON 


5 


• '/ -> o 


• u "0 


.07 0 


.071 


.0 80 


TARGET 


1 


WEAPON 


6 


.05 0 


. 0 0 


.100 


. 120 


. 1 30 


TARGET 


I 


A E A PON 


7 


.120 


. 1 60 


. 130 


.200 


.202 


TARGET 


2 


WEAPON 


1 


. 0 1 0 


.0 20 


. 02 4 


.030 


.035 


TARGET 


-> 

«r- 


WEAPON 


2 


• »< 2 C 


. 030 


.04 0 


. 050 


.060 


TARGET 


'/ 


A E A PON 


3 


.012 


.021 


. 03 0 


.0 36 


.942 


TARGET 


2 


A E A PON 


4 


.016 


.032 


.04 7 


.061 


.076 


TARGET 


2 


WEAPON 


5 


.02 2 


.039 


.053 


.065 


.07 5 


TARGET 




WEAPON 


6 


.00 7 


.012 


.01 6 


.019 


.921 


TARGET 




WEAPON 


7 


.050 


. 1 50 


.15 5 


.200 


.240 


TARGET 




WEAPON 


1 


. 050 


.0 80 


. 140 


. 1 70 


. 1 90 


TARGET 


3 


WEAPON 


2 


.094 


. 143 


.172 


. 196 


. 20 5 


TARGET 


Zj 


WEAPON 


3 


.06 0 


.090 


.15 0 


. 190 


.22 0 


TARGET 


5 


WEAPON 


4 


.036 


.0 70 


. 102 


. 131 


.161 


TARGET 


5 


W EAP0N 


5 


.119 


. 175 


.22 0 


.233 


.245 


TARGET 


3 


WEAPON 


6 


.970 


. 100 


. 160 


.200 


.2 30 


TARGET 
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WEAPON 
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.03 0 


. 120 


.130 


.220 


.250 


TARGET 




W EA PON 
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.0 1 1 


.020 


.028 


.035 


.042 


TARGET 


4 


W EAPGN 


2 


.019 


.036 


.048 


.056 


.064 


TARGET 


4 


WEAPON 


3 


.019 


.0 32 


.042 


.05 1 


.057 


TARGET 


4 


WEAPON 


4 


.00 1 


. 001 


.00 1 


.001 


.002 


TARGET 
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W EA PON 
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. 0 i 2 


.0 54 


. 06 L 


.082 


.090 


TARGET 
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WEAPON 


6 


.00 9 


.016 


.02 3 


.028 


.033 


TARGET 


4 


WEAPON 


7 


.029 


.057 


.066 


.080 


.090 



TABLE 3 

KILL PROBABILITIES 



% 
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These weapons' weights and aircraft payload and rack capacities 
result in approximately 400 possible loadings which must be considered. 
From amongst these 400, Table 4 lists the five loads with the greatest 
expected returns and for comparison the loading yielding the least. 

A loading designated 1200002 denotes one weapon of type one, two 



: type two, 


zero weapons 


of types three through 


: type seven 


. + 




LOAD NO. 


LOADING 


EXPECTED RETURN 


1 


0300002 


1.29514 


2 


1200002 


1.25048 


3 


0100102 


1.22353 


4 


1000102 


1.20670 


5 


3000002 


1.20409 


381 


2011010 


.38069 



TABLE 4 

LOADINGS AND EXPECTED RETURNS, 
DYNAMIC PROGRAMMING METHOD 



The expected returns from the optimal load are over three times 
higher than those of the worst. In fact the twentieth best load, 

0020300, had expected returns of only .77700. Thus the optimal loading 
yields returns more than 65% better than a loading which is better than 
95% of the possible loads. 

The above inputs and constraints have also been used in the sub- 
optimal model to generate a "best load." Weapon 7 is the best single 
weapon against targets 1 and 2, and weapon 5 is best against targets 
3 and 4. Thus the only loads considered in this formulation were 
0000500, 0000401, and 0000102. Table 5 lists six possible load con- 
siderations (including stick options) and the expected returns from each. 
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Loading 



E 



0000102 (Weapon 7 used as a stick of 2) 



.5766 



0000401 (Weapon 5 used as a stick of 4) 



.4458 



0000102 (Weapon 7 used as 2 sticks of 1) 



.3824 



0000500 (Weapon 5 used as a stick of 5) 



.3740 



0000401 (Weapon 5 used as 4 sticks of 1) 



.3320 



0000500 (Weapon 5 used as 5 sticks of 1) 



.1432 



TABLE 5 



POSSIBLE LOADS AND EXPECTED 
RETURNS, SUB -OPTIMAL METHOD 



The "best load" using the sub-optimal hand method is 0000102 
considering weapon 7 as one stick of two. Although only the principal 
weapon for each target lias been considered, note that weight constraint 
still allows the loading of one more of the lighter, less effective 
weapons (weapon 2 rather than weapon 1 since weapon 2 is dominant over 
the entire target spectrum) . The load thus generated corresponds to 
the third most valuable dynamic programming load, and, if used in 
accordance with the dynamic programming firing policy, it would be 
expected to yield a return of greater than 907 o of that given by the 
actual optimum load. 

With the best loading as generated by either method, the optimal 
strategy can be listed using the "POLICY" program of Appendix II. The 
results of this determination for loading 0300002 are shown in Figure 
1 along with a format for displaying the strategy compactly. The 
time axis is mission time remaining and entries such 2x7 denote that 
two weapons of type seven should be expended. 

As an example, as the mission begins, the subloading is 0300002. 

If the first target which is met is a type one, and is encountered with 
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SUBLOADING 0300002 



1x7 



1x2 



NO FIRE 



2x7 



3x2 



3x2 



— -j- 

.15 



2x2 



.05 



.10 



.20 



.25 



su 


fBLOADING 


0300001 


1x7 


NO FIRE 


1x7 


3x2 


3x2 


2x2 



.05 



.10 



.15 



.20 



.25 



<M 

X 


NO FIRE 


3x2 


3x2 


3x2 


2x2 


1x2 


NO FIRE 


1 1 | 





SUBLOADING 0200002 



1x7 



1x2 



NO FIRE 



2x7 



2x2 



ix7 



2x2 



.05 



.10 



.15 



• 20 



25 



FIGURE 1 

OPTIMAL STRATEGY DISPLAY 
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TARGET 
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3 



o 
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SUBLOADING 0100000 
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SUBLOADING 0000002 
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2x7 
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.25 



1x7 


NO FIRE 


1x7 


1x7 


1x7 

r 


NO FIRE 

i l I 



.25 



(The abscissa in these figures is the time remaining in which to 
search for and attack targets.) 
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.22 hours remaining, the optimal strategy indicates that no weapons 
should be fired. If the next target which is encountered is a type 
four at time .17, the subloading still being 0300002, a drop of three 
weapons of type two is called for. The sub-loading now is 0000002. 

If a type three target is now met, both of the remaining type seven 
weapons should be dropped. 

It is of interest to note the changes which occur over time in 
the optimal weapon to drop when a target is met. With the loading 
0300002 and target type one, with more than .08 hours remaining, no 
weapons should be expended. With .045 to .08 hours remaining, the 
optimal strategy calls for the expenditure of one weapon of type two. 
Weapon type seven however has a higher kill probability against type 
one targets. 1 The type seven weapons at this point are being saved 
for possible future encounters with targets of type two and four, both 
of which are more valuable and occur frequently enough so that there 
is a significant probability of an encounter in the time remaining. 
Although weapon type two is the most potent against type three targets 
-the most valuable target type- the probability of meeting a type 
three target is rather small because of the low encounter rate of 0.10. 
It is for this reason that type two weapons are not being conserved. 
Later in the mission however, the optimal weapons to drop on a type 
one target are type sevens. At this point in time the probability of 
future encounters with more valuable targets is remote. 

Since the ultimate objective is the maximization of the expected 
returns for the mission, a consideration might be made for lengthening 
the search time at the expense of lower payload and weapon station 
availability. That is, perhaps an additional tank of fuel can be 
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carried on one of the aircraft's weapon racks, thus reducing the 
number available on which ordnance can be carried and also reducing 



the weapons payload, but increasing the search time* 

Table 6 lists the expected returns for the five best loads 
generated when the search time has been increased to 0.50 hours, the 
payload decreased to 4000 lbs., the number of weapon racks decreased 
to four, and the other input parameters unchanged from those of Table 
3. 



Since the expected value of targets killed for the longer mission 
is significantly better than that for the shorter, 1.78803 compared to 
1.29514, it is clearly preferrable to make the trade-off of gas for 
weapons . 

An interesting feature of the kill probability matrix of Table 
3 is the preferred weapons with which to attack the most valuable 
targets. Weapon type seven is the most potent against type two tar- 
gets while weapon type five is preferred against type three targets. 
Target type three, although the most valuable, has a very low encounter 
rate. This results in the fact that weapon type five does not appear 
frequently in the loadings which yield high expected returns for short 
missions. Letting the maximum payload equal 4000 lbs. and the number 
of weapons stations equal four, the preferred loadings of these two 



LOADING 



EXPECTED RETURNS 



0000002 

0000301 

0100201 

0013231 

0000400 



1.78803 

1.42184 

1.30025 

1.29743 

1.28808 



TABLE 6 

BEST LOADS AND EXPECTED RETURNS 
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weapons alone can be examined as a function of search time. The three 
loadings which must be considered are 02, 31, and 40, where now, 02 
denotes a load consisting of zero weapons of type five and two weapons 
of type seven. Calculations show that an initial loading of 02 is 
optimal for mission times of less than 1.6 hours. For mission times 
of 1.6 to 8.9 hours, loading 31 is optimal, and finally for T > 8.9 
hours, loading 40 is preferred. 

COMPARISON OF EXPECTED RETURNS FROM NON-OPTIMAL 
FIRING STRATEGIES 

The gains in kill potential realizable through the utilization 
of the optimal firing doctrine can be demonstrated by comparisons with 
the expected returns from other strategies. As an example, consider 
the following case with non-optima 1 strategies: 

Let: 

Number of weapons * 1 
Operating time *= 1 hour 
No. of target types = 2 




v =4 
2 

Strategy A: Expend the weapon on the first target which is 

encountered . 

E = the expected returns from strategy A. 

A 

Strategy B: Expend weapon only if a target of type 2 is 

encountered. (Note that target type 2 is the 
more valuable). 

E = the expected returns from strategy B. 

B 

T.: The random variable describing the time to first encounter 

1 with target type i. (T . is exponentially distributed with 
mean 1/r^.) 1 
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E = p v.Prob. (T < min (t , 1 hour)) + p v Prob. (T < min (t , 1 hour)) 
A* 1 1 2 22 2 1 



1 -r t “ -r t 1 -r t 

= P n v re lire 2 2 dt dt + p 0 v re 22 

^ 1 1 J l «J 2 21 2 2 J 2 



-r t 

r e 1 1 dt dt 

1 1 2 
t 2 

E - p v Prbb. (T < 1 hour) 

B 2 2 2 

1 

-r t 

K = P 9 V 9 I r o e 22 dt 
B 2 2 J 2 2 

o 

For several values of the parameters p. and r , Table 7 shows the 

1 i 

results of calculations of the expected returns using either of these 
strategies or the optimal. 
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TABLE 7 

GAINS IN EXPECTED VALUE OF TARGETS KILLED ACHIEVED BY 
UTILIZING THE OPTIMAL FIRING STRATEGY 
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These results are representative of those obtainable using other 
values for the input parameters. Generally, when one of these non- 
optimal strategies approaches th 2 optimal in expected returns, the 
other yields very poor results. When the non-optimal strategies yield 
roughly equivalent returns, the optimal dominates both by a significant 
amount; this is demonstrated in the first example of the table by gains 
of 267. and 327* through use of the optimal strategy. 

Comparisons can be made in this way with any conceivable strategy. 
It seems reasonable too that techniques which consider the target 
values and encounter rates could yield returns comparable to those 
achievable through use of the optimal strategy. The determination of 
these strategies however would probably involve calculations no less 
complex than those of the dynamic program. 



SENSITIVITY OF EXPECTED RETURNS TO TARGET VALUES 

A complete parametric analysis is beyond the scope of this paper; 
however, some further comment must be made about target values. The 
developments of both the dynamic programming model and the hand appro- 
ximation are based upon the assumption that a ratio scale of values can 
be assigned to all possible targets. Such a determination is completely 
subjective and must rely solely on the judgement of experienced planners. 
Quantitative target values might be determined by presenting a list of 
targets to several military commanders and asking them to rate these 
targets according to military value. Targets could be assigned numbers 
consistent with their importance relative to a target of least value. 
The results of such a large number of opinions could then be analyzed 
by methods of sensory psychophysics in order to assign representative 
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values to all targets* Methods for accomplishing such an analysis 

have been reviewed and consolidated by S. S. Stevens.'*' 

The above method will obviously generate strategic values. These 

will certainly not be valid for every tactical situation (e.g., a tank 

about to overrun friendly troops is more important than a railroad 

bridge at that moment even though the bridge may have a higher strategic 

value). A strategic scale will, however, provide a starting point -- 

a base which can be adjusted to meet specific tactical situations. 

An analysis of the effects on the expected returns of incorrect 

selections of target values can be made rather easily for the simple 

case of a single weapon and two target types. Assume the "true' 1 ratio 

of the target values is V *• Now suppose that an erroneous ratio, 

V r , has been used to calculate the strategy. The expected returns 

which result from the utilization of this strategy - these returns 

being computed using the true value ratio, V * - can be compared with 

r 

the maximum of the expected returns, this maximum being obtained when 

the strategy derived using V * is utilized. Define: 

r 



V r = V 2 /Vj 

and let the mission time - 1 hour. With a single weapon, a strategy 
consists merely of a cross-over time, the time after which the least 
valuable target should be attacked if met. These times can be calcu- 
lated from: 

r lPl + r 2 p 2 v r 

t = — L_ i n Z r if V > 1 



r l +r 2 r 2 (p 2 V r - p l> 



t = 



1 r ' 



In 



W + r 2 p 2 V r 
r l (p l - p 2 V r > 



t = 00 



if V < 1 
r 



if V =1 
r 



o.S. Stevens, "A Metric for Social Concens us , V Science , vol. 151, 
4 February 1966, pp. 530-541. 
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Using these cross-over times then, the expected returns, E, are cal- 
culated from: 



r 9 -< r l +r 2 ) P i r l ‘ (r l +r 2 ) 

E - p^V_* — -L- (1 - e ) + riT - (1 - e z )lftsl hour 



2 r r + r 
1 2 



r +r 
1 2 



-r (1-t) -r 2 (l-t) -(r +r )t r 

E “ p V * (1 - e L + e ( 1 - e ) — — ) 



2 r 



r i +r 2 



+ p e 2 (l- e l2 ) if t < 1 hour, and V > 1 

*1 v 7 r +r r 

1 2 



E - p i (l - e 



" r i ( 1-t) -r,(l-t) -(r fr,)t 



+ e 



12 

(1 ' e 



-r (1-t) -(r +r )t r 

+ p V * e 1 (l-e 1 2 ) — - — if t < 1 hour, and V < 1 

2 r v } r^+r 2 r 



Letting: 

p “ p = . 10 

1 2 

r “ 2/ hour 
1 

r = 1/hour 
2 

Table 8 lists the results for several possible values of V *. 

r 

While the scope of these investigations of value sensitivity is 
obviously limited, it appears that over-estimation of the value of 
targets of higher worth might impose less risk than under-estimation. 
Further examination in the area of sensitivity where several weapon 
types are employed simultaneously would be necessary before conclusive 
results could be stated. 
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If V * is misestimated by a factor of z, 




r 










the expected returns decrease by y. 
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<.57. 
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TABLE 8 

PERCENTAGE DECREASE IN EXPECTED RETURNS 
DUE TO MISESTIMATE OF V * 
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4. CONCLUSIONS AND RECOMMENDATIONS 



A method for determining an optimal initial loading for missions 
against targets of opportunity has been developed. As a by-product 
of this load determination, an optimal strategy for the sequential 
expenditure of weapons has been found. The gains in the expected 
value of targets killed, when the initial loading is optimal and when 
an optimal strategy is employed, have been shown to be significant 
when compared with initial loadings and firing doctrines derived by 
the application of simple intuition and judgement to kill probabilities, 
target values, and encounter rates. 

The determination of an optimal loading and expenditure policy, 
however, are at present of such complexity as to require the use of a 
high speed computer in all but the most trivial cases. For those 
operational situations in which a computer will not be available, a 
hand-calculable approximation method for the determination of a best 
initial loading has also been developed, which, though sub-optimal, 
has been found to generate loads which yield high expected returns. 

Either of the above methods should yield significantly better 
results than present policies which prescribe that aircraft be loaded 
for targets which are assumed known or that several aircraft each be 
given different loads in order to cover an unknown spectrum of targets. 

The models presented in this paper assume that weapons can be 
loaded in any combination which satisfies rack and weight constraints, 
and that the pilot has complete freedom of choice as to the sequence 
in which weapons can be expended. For most attack aircraft these 
assumptions are not always true. Restrictions on unsymmetrical loadings 
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are one limiting factor. Others relate to the mechanics of multiple 
bomb racks. The models should be expanded to account for these aspects 
of the problem. 

A vital assumption in the developments of both models was that 
inter-encounter times would be distributed exponentially. It is 
recommended that an analysis of operational data be made to ascertain 
the validity of this supposition. 

It seems reasonable to assume that a list of all possib3.e types 
of targets of opportunity, for a spectrum of scenarios, would be 
small enough - certainly not larger than 50 - that a thorough analysis, 
as suggested in Section 3, could be made of their relative values. Such 
an analysis would indicate general strategic values which could then be 
adjusted for specific tactical scenarios. Furthermore, additional 
investigations similar to that outlined in Section 3 should be under- 
taken to study the sensitivity of the optimal load and strategy cal- 
culations to these values. 

The hand approximation method is rather crude in form. It is 
recommended that this model be examined further in an attempt to refine 
the calculations and to extend the model to cover a wider variety of 
cases in a more complete fashion. 

The computer programs of Appendix II were written solely for the 
computations used in this thesis and should be improved upon before 
extensive investigations are made of the dynamic programming treatment 
of the problem. The four target-seven weapon example required two 
hours and fifteen minutes on a CDC-1604 computer for the determination 
of an optimal load. This is testimony to the inefficiency of the 
"BESTLD" FORTRAN program. 
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The problem of determining an optimal initial loading is similar 
to the classic "flyaway-kit 11 problem 1 * of dynamic programming. In the 
strike aircraft loading problem the additional complication is present 
that the costs of shortages are not constant but depend on the initial 
loading. It would seem though that research into solution of the 
strike aircraft problem by this method should be fruitful. Very good 
approximations to an optimal initial loading should be obtainable very 
rapidly in this way. 

Although the "POLICY" FORTRAN program computes the optimal 
strategy with relative ease, usually less than four minutes is required 
on a CDC-1604 computer, it is recommended that methods of approximation 
be developed for which a strategy can be calculated without the aid of 
a high speed computer. 



1 

Bellman, R. E. and Dreyfus, R. E., Applied Dynamic Programming , 
(Princeton, New Jersey: Princeton University Press, 1962), pp* 42-43. 
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APPENDIX I 



ENCOUNTER RATES 



In a combat situation, many targets (especially "targets of 
opportunity 11 ) may be thought of as being distributed randomly through 
out a given area, in that the pilot flying the mission does not know, 
beforehand, where or when a target will be encountered. This may be 
thought of as the classical situation in which it may be assumed that 
in any small time interval dt, the probability that a target is en- 
countered is proportional to dt and is independent of any previous 
encounters. 

Over a period of time, photo or visual reconnaissance can give 
an average density of targets in a given area. From the attack air- 
craft speed and the average search sweep width (width of a path in 
which the pilot could be expected to see a target) the area swept per 
unit time can be obtained. By multiplying the average type i target 
density by the average sweep area per unit time, an estimate of the 
average encounter rate r^ is obtained for target type i. 

The average encounter rate developed above can be used to gen- 
erate the distribution of inter-encounter times. 

Assumptions: 

Probability of one encounter in dt - r^dt + h(dt) where h(dt) 
approaches zero as dt approaches zero. 

Probability of two or more encounters in dt = h(dt) 

Probability of no encounters in dt = 1-r.^dt + h(dt) p(t)dt = 

prob. (t ^ first encounter time ^ t + dt) = dP(t) 

rt 



P(t) 



p(T)dt = prob. (encounter time ^ t) 



o 
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Q(t) = l-P(t) 



Thus if t is the encounter time, and t is any arbitrary time in 
the mission, 

Q(t+dt) = prob.((t > t)(no encounter in dt)) 

e 

- Q(t) (1-r.dt + h(dt) ) 

Subtracting Q(t) from both sides, dividing through by dt and taking 

the limit as dt approaches zero gives dQ(t)/dt = -Q(T)r 

Rearranging gives dQ(t)/Q(t) = -r dt 

i 

Integrate both sides to get In Q(t) = C-r t 
Thus Q(t) ■ Ae r i t 

But Q(0) = 1 since there can certainly be no encounters in zero search 

-r.t -r.t 

time. Therefore A * 1 and Q(t) = e 1 and consequently P(t) =*l-e 1 

— r r 

and p(t) = r^e i 

Thus it is seen that the inter-encounter times are exponentially 
distributed. It also follows naturally that the probability distribution 
on the number of targets type i encountered during a mission of length 
T is Poisson, ^ i.e. , 

n -r .T 

prob. (N(T) * n) * (r^T) e 1 /n! 

and the expected no. N(T) = r T (This is obtained from the observations). 



Parzen, E., Modern Probability Theory and Its Applications , 
(New York: John Wiley and Sons, 1960) p. 253. 
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BEST LD AND POLICY FORTRAN PROGRAMS 
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XLAM(I) IS THE ENCOUNTER R A TF WITH TARGET TYPE I* 000180 

V ( I ) IS THE VALUE OF TARGET TYPE I 000190 

DO 825 K= 1 ♦ M 000200 

DO 825 JN = 1 ♦ LL 000210 

READ 826, ( P ( I , JN , K ) , I = 1 , L ) 000220 



o 


o 


o 


o 


o 


O 


O 


o c 


o 


o 


c 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


c 


o 


o 


o 


o 


o 


o 


o 


o 


o 


CO 




in 


vo 


r~ 


c c 


a 


o 




(NJ 


r< 




lT\ 


v£5 


r- 


CO 


a 


o 


f — < 


CM 


CO 


<r 


in vO 


r~ 


00 


o 


o 


r— » 


r\j 


CO 




CM 


<\J 


cm 


r\i 


r\j 


cni 


(\j 


cr 


co 


co 




CO 


co 


co 


co 


CO 


CO 


<r 


<r 




<* 


<r 






<r 






m 


m 


in 


uo 


n 


O 


o 


o 


o 


O 


O 


O 


o 


o 


O 


C 


c 


O 


O 


o 


O 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


a 


o 


o 


o 


o 


o 


O 


Q 


O 


o 


c 


O 


O 


o 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


Q 


O 


o 


c 


c 


c> 


o 


o 


a 


o 


o 


o 


o 


o 


c 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 





























«—• 
































































X 












II 




























< 








9 
















X 
















































X 
















CO 












m 




























x 








vO 






X 










• 












e 




























x 








• 






X 










in 












0 




























•— 








n- 






X 










X 












X 




































X 






\ 










» 












< 




































*-H 






• 










in 












X 




























_J 








•• 






C 










•— i 












> 




























J 








X 






• 














— 




































•— 








vO 






vO 










y— 




X 








X 




































9 






Ll 














X 








0 




































vO 






* 










< 




U 




































z 








X 






CO 










y— 




< 








X 




























“) 








»— 4 






*— • 










_J 




> 








X 




















































X 




X 








CO 




























X 








X 






e 










O 




in 




































X 








v£> 






< 














#- 








X) 




























>- 








* 






o 










X 




X 




































»— 








lT 




























X 




































X 






X 










*- H 












0 










— 


















X 








»—♦ 






<r 








X 






X 








r— 4 










CO 


















O 








• 






a. 








♦ 


X 




1 — 


















» 


























X 


— 












(NJ 


*■ 




<c 








X 










vO 


















in 








v£> 


\ 




> 










in 




X 








vO 










X 


















z 








• 


\ 




< 










X 






















co 


















o 








<r 


\ 




y 








in 


x> 




X 








N 










*■ 


















a 








X 


X 














o 




LU 








#■ 










X 


















< 








. — i 


c 




X 










X 




H- 








in 










N, 


















UJ 








*• 


• 




lT' 








in 














X 










#■ 


















:< 








X 


r- 




f — * 








UJ 


X 




0 








•— > 










(NJ 


























o 


Ll 




» 








•a 


vO 




0 








h- 










4- -4 


















•— i 








*• 


00 




(Nj 








X 






u 








4 — . 










»■ 


























CO 


» 




•— 








h- 


(Nj 




z: 


— 






X 








4— 


z 


















►— 








X 


X 




* 










• 




X 








*— 








X 


0 


















< 








•-H 


L 




cO 








x 


in 




X 








cr 










X 


















X 








r 


X 




*— > 








UJ 


X 






*— t 


— 




< 








r— 1 


<r 


















f— 








X 


O 












O 


*■ 




r-“ 4 


II 


(NJ 




CC 








ll 


X 


























vO 


»--• 




in 








IX 


c0 






— * 


• 




0 








•— 




















> 








*4 


UJ 




v 








< 


•— • 




X 




iT> 




X 










X 


















h- 








(NJ 


y 




N' 








♦ — 






<r 


— 


X 




X 








— 


vO 


















— > 








X 






< 










X 






»— 


» 














#> 


















x 








f—i 


ZZ 




CL 








Ll 


V 




X 




X 




X 








*• 


X 


















— ■ 






— * 




o 












O 


— ■ 




X 


> 


vO 




X 








z 


CO 


















QD 






00 


X 


Cl 




U_ 










h- 




X 


+ 


+ 




*— i 








~) 


#4 


















< 






#4 


vO 


< 




O 








cr. 






1 — 


— 


lT 




V 








#> 


CO 


















GO 






f-H 


♦ 


Lu 










< 


UJ 


2Z 






►-4 


» 


— 


X 








*— 4 


•—4 


















o 






II 


I 


5 




CL 






y— 


CC 


O 




t— 




c 


X 


00 


— 








•- 


















< X 






1 — t 


X 






Ll 






_J 


V 






X 


ZF 


•— l 




r— l 


vO 






X 


h- 










— « 








x 


• 




*4 


f-H 


X 




cr 






LU 


b 


in 




0 


< 


X 


H 


r 










X 










4—4 




O 






V 




— 




in 




>' 






C 




in 




X 


X 


#- 


II 


X 


r- 






#> 


0 
















— 


UJ 








X 


f— 4 




3 








X 


*— • 




< 


X 


X 


•— 


r-* 


+■ 




x 2 


X 














r-H 


o 


X 


UJ 




— 


CO 


X 


Z\ 


X 






h- 


in 






y— 




c 


#> 


(NJ 


r-H 




X 


“0 


< 










♦- 


0 




• 


y— 


a 






f— 4 


* 


< 










(NJ 








— 


f~4 


*— 


*4 


•— 4 




#4 


#4 


y— 






— 


co 


— * 


• 


II 


in 




> 








oo 


9- 


X 








e 


x 












N 


#4 


#> 


f— 4 










»-H 


#4 




» — 1 




Ll 


ii 


y— 




#4 


f— “4 


X 


X 


O'' 








X 


vO 




(NJ 


*4 






X 


X 


rH 


ll 




X 






X 


(NJ 




+ 


X 


00 








CO 


X 


*— H 










<r 


r-H 


f— * 


in 


* * 


vO 


►— . 


r- 


\ 


<r 


II 


2: 


vO 


I s - 




CO 


f-H 


II 




< 


f"H 


> — ■ 


— 


X 


UJ 


in 


r— < 




in 








o 






0 




0 


— 


c^ 


— 






“> 


CO 


— 


X 


CO 


w 


*—4 




1 


”0 






O 


r> 


in 






in 








CO 






CO 




CO 




00 










00 




X 


r- 






X 






»— 






x 




h- 




in 


h- 


o 






h- 






h- 




»— 




h- 




in 


in 




h- 


z 




h- 


O 


1 — 1 


> 


i-H 


<• 


z 


i— 




y— 


<r 




I— 


< 


• 




X 


< 




K— 


< 


y— 


< 


i— 


< 




CO 




y— 


< 


4— 4 


h- 


< 


O 


z 


11 


O 




■o 


UJ 


h~ 


X 


X 




Z 


V 


o 


O 


X 






X 




z 










00 


00 Z 


21 


1— 


z 


"C - 


n 


1 4 


z 


in 


X 


#> 


o 


Z 


*— ■ 


(X 




*— < 




m 


II 




X 




*—4 


X 


• — < 


X 


♦— < 


X 








> — « 


X 


z 


•— i 


X 




IF 


•—4 




o 




x 


o 


X 


o 




X O 


X' 


in 


cr 


o 




X 


0 


X 


O 


X 


0 




O 


O X 


0 


0 


cr 


0 


O 


ll 




C 


li- 




< 


u 


CL 


Ll 




Cl 


Ll 


M 




CL 


X 




a 


X 


X 


X 


X 


X 




O 


C X 


X 


U X 


X 


O 


< 


< 


O 




CL 


i— 








H 














i-H 














r-4 
























ve 






in 




CO 






m 














in 




vO 




X- 










vO 


in 




co 




O 






nj 






(NJ 




in 






in 








o 






O 




O 




co 










CO 


CO 




CO 




O 






00 






00 




in 






in 








CO 






CO 




CO 




00 










00 


00 




n* 




in 







u u 



53 



ooooooooooooooooooooooooooo 
m vd r- co o o *— ir\jco<fo^vDf^ooc^o*— if\iro<}-irio( s -ooo , 'OH 

ir\ir\ir\inin\D'0v0'0'Dv0viD\DviDv0r s -( s -(^i > *f > ‘( s *f^r^(^f s “COaD 

ooooopoooooooooooooooooooco 

ooooooooooooooooooooooooooo 

oooooOooooooooooooooooooooo 











F- 






F- 






1 — 




F- 




1 — 






F- 












c\j 






co 












in 




O 




I s * - 






00 


















“> 






“> 






> 




“> 




"> 






“> 












#s 






•> 






e- 






e- 




*> 




* 






•n 












rH 




— 


rH 




— 


rH 






rH 




— * rH 




*—* rH 














00 




II 




f\j 


II 




CO 


II 






II 




m 11 




vD II 




r- 


11 










H 


— CL 


rH 




CL 


r-H 




CL 


*— 1 


— 


a. 


*-h 


— Cl 


rH 


— CL 


rH 




a 


«— 1 






rH 


i 


rH C\J 


1 




CO 


1 


“> 




1 


“O in 


I 


“> O 


1 


“> r- 


1 


”> 


00 


1 


0 


0 


II 


CL 


— CL 


l 


”> CL 


1 




CL 


1 


“> Cl 


1 “> 


CL 


1 -> 


CL 


1 


■> 


CL 


• 




k— < 


rH 


"> 


CM 


F- 




CO 


1 — 






i— 




in 


F- 


vD 


i— 


n- 


F- 




00 


c 


II 






I 




cm 


rH 


“> 


CO 


rH 






«— » 




in * — 1 


“> 


vO H 




r- 


rH 


-> 


II 




nj 


11 


O' O 


II 


~> 


O 


11 


“) 


O 


II 




O 


11 


“> 0 


II 


“> O 


11 


“> 


O 


II 


0 


0 


O 




11 in 




11 


m 




II 


in 


— ■ - 


11 


in 




11 in 


«— 


11 m 




11 


m 


— * 


< 


< 


m 


t—i 


h- 


cm 


h- 




CO 


F- 






1— 




m 


h- 


O 


I— 


n- 


1— 




00 


0 0 






nj 0 




co 


O 


■ — 




O 




m 0 


■ — 


0 O 




I s - O 


w 


00 


C 


■ 


1 


_J 


O 




•> O 






0 


-> 


~) 


Q 






O 




"> O 




“> Q 






O 




3: 


2: 


O 



— ^ 
* 

4- •— • 

o -> 
< + 
c o 

-J < 

2: o 

-J 
II 3: 
II 

o o 
< < 
o o 

-J -J 

2: 5 



r\J 

o 

in 



54 



IF (NLOAD «LE* L) 823*501 000820 

823 IF { WLOAD .GE. AMIN) 503*506 000830 

506 IF (NLOAD »EQ* L) 504*501 000840 

503 IF(WL0AD »LE. WM) 504*501 000850 

504 NIP = J(l) + 1 000860 
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18 = NN8-1 001140 

JPROD =11 * 10*»7 +12 * 10**6+ 13 * 10**5 + 14 * 10**4 001 150 

1+ 15 * 10**3 + 16 * 10**2 + 17 * 10 + 18 001160 

JJ = JJ + 1 001170 

1000 ISUBS(JJ) = JPROD 001180 



LOAD = ISUBS(IPROD) 001190 

58 SUMLAM = 0.0 001200 

DO 100 I = 1 » M 001210 

100 SUML AM= SUML AM + XLAM(I) 001220 
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DO 108 16 = 1 * I PROD 001510 

IF ( JSUBS - I SUBS (16)) 108,109,108 001520 

109 17 = 16 001530 

GO TO 110 001540 

108 CONTINUE 001550 
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FOOOOO 1 20 ( .050) = .34921 
FOOOOO 1 2 0 ( . 100) = .56935 
FOOOOO 1 20 ( .150) = .76357 
FOOOOO 1 2 0 ( .200) = .94750 
F00000120 ( .250) = 1.12276 
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PROGRAM POLICY 000000 

DIMENSION XLAM( 10) * P ( 8 * 8 ♦ 1 0 ) * V ( 10) ♦ X ( 5 1 2 ) ♦ F ( 512) > N ( 8)* 000010 

1 I SUBS ( 512) »KFIRE(9»lO) * JF I R E ( 9 ) * L F I RE ( 9 ) 000020 

RFAD 200»M,N1*N2*N^,N4*N5»N6»N7»N3.DFLTA,T 000030 

200 FORMAT ( 9I3.2F5.0) 000040 
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N5P=NS+1 000270 
N6P=N6+1 000280 
N7P=N7+ 1 000290 
NRP = N8 + 1 000300 
IPROD = NIP * N2P * N3P * N4P * N5P * N6P * N7P * N8P 000310 



DO 101 I=1»IPR0D 000320 
X { I ) = 0.0 000330 
JJ=0 000340 
DO 1000 NN1 = 1 » N 1 P 000350 
II = NN1 - 1 000360 
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NL = N ( 14 ) + 1 000960 

IF (NL-1) 105*105.107 000970 

DO 106 I 5 = 1 * NL 000980 

J I 4 =15-1 000990 

J 1 4 IS THE NUMBER OF WEAPONS OF TYPE 14 EXPENDED 001000 



oc oooooooooooooooooooo 
HOsjc^vj-in^r-coo o*-H<\jro<j'cn\0( s ^coa^OrHfvj 

OOOOOOOOOrHrHrHrHrHrHrHrHrH,— lOsjCMCM 



oooooooooooooooooooooo 

oooooooooooooooooooooo 















CO 






























< 














o 






























►— 














rH 






























-J 












































LU 












































O 










*— « 




O 






























* 










1 




rH 






























r— 








CO 


CO 




» 






























CO 








*— t 






00 






























•— < 








w 






o 






























— 








> 






•H 






























5 : 






































CO 






< 








* 


o 












































rH 


































X 










* 




v£> 




















CM 












rH 






CO 




O 


•— i 












rH 








rH 




t — * 






* 








*— I 


*— i 


O 


■w 












rH 








rH 




— » 








vO 








3 


oc 


if) 












rH 




3 




» 




LU 








♦> 






vf 




CL 


CD 










O 










CM 




CL 








o 






*— t 


1 




3 










z: 


rH 




It 




rH 




>— 1 






+ 


<1- 










* 


to 










0 


rH 








rH 




Ll 








vD 






<}• 


o 


rH 


•— « 










u 


rH 




— - 




* 




X 






CM 








•— < 


< 












r— 


LU 






CO 




in 










51 


O 


o 




3 


o 


11 


1 








I s - 


to 


vO 




»— 1 




0 




ll 


*-* 




CL 




• 




*— ■ 


-J 












*— 1 




O 








rH 






11 




LU 


vO 


o 


CM 


CL 




vO 


to 




o 






+ 


rH 












— * 






I'- 












*— 1 


CD 




rH 


LU 


X 








♦— * 


LU 


— ' 




CO 


CO 


LU 






II 




11 


11 




3 


'O 


rH 


3 


11 


f— 


< 




— 


3 


>- 




•— t 


*— 1 


3 


ll 












CO 


to 


t— » 




2 


o 


to 


1 






Z 


1 






• — 


e 




*— < 


I- o 


1— 


to 


o 


3 




o 


M 1 


2 : 


CL 


rsl 


M 


LU 




< 


< LU 


LU 


*— < 


CM 


3 


to 


1 - to 


CD 


rH 


• — 


11 


1— 


1— 


0 


*— < 






CL 


t— 






CL 


CL 


1— 


£ 


— 


cl 




CL 


3 










e 


u 


Ll 




It 


*— « 






If 


*— « 


*—t 




CL 


Ll 


t—* 


o 




to 


O Li- 


r- 


o 


o 


LU 


11 


Ll 




Ll 


O 


Ll 




Ll 


Ll 


O 


LU 


•— < 


Ll 


C Ll 


“O 


O 


•— 1 


_• 


o 


u 


CO 


M 




< 




u 




X 3 




U 


1— 




O 




i — i 


C\l 






O' 




00 


c 






rH 




vO 




CM 






iTl 






vj- 












o 




o 


rH 






rH 




0 




rH 






0 


O 




vO 




vO 


vD 






rH 




rH 


rH 






rH 




H 




rH 






rH 


»H 



6 2 



F ( I 2 ) = T ERM 1 + TERM2 001230 

PRINT 302.(N(I) . 1 = 1,8) 001240 

302 FORMAT <15H SUBLOADING IS , 8 1 1 ) 001250 

PRINT 303, ( I » JF I R E ( I) » LF I RE ( I ) » I = 1 »M ) 001260 

303 FORMAT (15H IF TARGET TYPF,I2,20H IS ENCOUNTERED FIRE.I2. 001270 



1 8H OF TYPE » ! 2 ) 001280 

103 CONTINUE 001290 

PRINT 320, ( N( IY) , I Y= 1 , 3 ) *T IME, C ( IPROD) 001300 

320 FORMAT ( / / , 1 X , 1 HF , 3 I 1 , 1 H ( , F 5 . 3 , 3H ) = , F 1 0 . 5 ,////// ) 001310 

DO 113 1=2, IPROD 001320 
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TABLE 10 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

BJ 
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2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 



APPENDIX III 



TABULATED VALUES OF 



* .05 


.10 


.15 


O 

CM 

• 


.25 


.04877 


.09516 


' .13929 


.18127 


.22120 


.04998 


.09904 


. 1 4948 


.19879 


.24770 


.05000 


.10000 


.14998 


. 19994 


.24986 


.05000 


.10000 


. 15000 


.20000 


.24999 


.05000 


. 10000 


. 1 5000 


.20000 


.25000 


.05000 


. 10000 


. 15 C 00 


.20000 


.25000 


.05000 


.10000 


.15000 1 


.20000 


.25000 


.05000 


. 10000 


. 15000 


.20000 


.25000 


.05000 


. 10000 


• .15000 


.20000 


j .25000 


.05000 


.10000 


. 1 5000 


.20000 


1 .25000 


.05000 


.10000 


.15000 


.20000 


.25000 


.05000 


. 10000 


.15000 


.20000 


.25000 


.05000 


. 10000 


.15000 


.20000 


.25000 


.05000 


.10000 


.15000 


.20000 


.25000 


.05000 


.10000 


. 15000 


.20000 


.25000 


.05000 


.10000 


.15000 


.20000 


.25000 


.05000 


.10000 


.15000 


. 20 G 00 


.25000 


.05000 


. 10000 


.15000 


.20000 


. 25000 . 


.05000 


. 10000 


.15000 


.20000 


.25000 


.05000 


.10000 


.15000 


.20000 


.25000 



RT 



= .30 

.25918 

.29612 

.29972 

.29998 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 

.30000 



.35 

.29531 

.34398 

. 34949 

.54997 

• 35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 

.35000 



.40 

.32968 
.39123 
. 3991 6 
.39993 
.40000 
.40000 
. 40000 
.40000 
.40000 
• 40000 
.40000 
.40000 
.40000 
.40000 
.40000 
.40000 
.40000 
.40000 
.40000 
.40000 



.45 



.43781 

.44869 

.44989 

.44999 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

. 45 C 00 

.45000 



.50 

.39347 
.48567 
.49806 
.4998 1 
.49998 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
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VALUES OF f . . 
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0 J 

1 

2 

3 

4 

5 

6 

7 

R 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

BJ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

1 3 

1 4 

15 

16 

17 

18 

19 

20 



= .55 


.60 


.65 


.70 


.75 


.42305 1 


.45119 


.47795 


.50341 


.52763 


.52878 


.57309 


.61658 


.65922 


. 70099 


.54724 


.59621 


. 64492 


.69336 


.74150 


. 314 97 1 


.59956 


.64937 


.6991 1 


.74879 


. 54997 


.59996 


. 64993 


.69990 


. 74985 


. 5 5 j 0 0 


.60000 


. 64999 


.69999 


.74998 


. 55000 


. 6000 0 


.65000 


. 70 ( 00 


.75000 


.55000 i 


.60000 


.65000 


. 70000 


.75000 


.55000 i 


.60000 


.65000 


.70000 


w 7 5000 


.55000 1 


. 60000 


.65000 


.70000 


4 75000 


• 5 5 0 0 < J 


. 60000 


.6 5 C 0 0 


.70000 


.75000 


. 55000 


.60000 


.65000 


.70000 


.75000 


. 55 J 00 


. 6000 u 


.65000 


. 70 C 00 


.75000 


. 55000 


.60000 


.65000 


.70000 


.75000 


. 55000 


. 60000 


.65000 


.70000 


.75000 


. 55000 


.60000 


.65000 


.70000 


. 7 5 C 00 


. S 5 uOO 


. 60000 


.65000 


.70000 


. 75 COO 


. 55000 


.60000 


.65000 


. 70000 


.75000 


. 55000 


.60000 


.65000 


.70000 


.75000 


. 55000 


.60000 


.65000 


. 70 C 00 


.75000 



RT = .Co .85 



. 5506 7 


. 57259 


.74188 


.78187 


. 78930 


.83675 


.79838 


.84788 


. 79979 


.8497 1 


.79998 


.84997 


. 80000 


.05000 


.80000 


.85000 


. 80000 


.85000 


. 80000 


.85000 


.80000 


.85000 


. 80000 


. 850 00 


. 80000 


.85000 


.80000 


.85000 


.80000 


. 8 5000 


.80000 


.85000 


.80000 


.85000 


.80000 


.85000 


.80000 


. 850 00 


.80000 


.85000 



. 9 -j 


.95 


.59343 


.61326 


.62095 


.6591 1 


.68380 


.93045 


.89726 


. 94653 


.89961 


.94948 


.89995 


. 94993 


. 89999 


.94999 


. 90 C 00 


.95000 


.90000 


. 950 JO 


' .90000 


.95000 


.90000 


. 9 5000 


. 90 G 00 


.95000 


.90000 


.95000 


.90000 


.9 5 ' 00 


.90000 


.95000 


. 90000 


.95000 


.90000 


.95000 


.90000 


.95000 


.90000 


. 95 C 00 


.90000 


.95000 
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BJS 


r r x i 


2 


1 


. 63212 


..86466 


2 


. 89636 


1.45866 


3 


, .97666 


1 .78198 


4 


. 99565 


1.92486 


5 


t * 99931 


1.97751 


6 


. 99991 


1 . 99408 


7 


. 99999 


1 . 99861 


6 


1.00 000 


1 . 99971 


9 


1.0 00 0 0 


1.99994 


10 


1.00000 


1.99999 


11 


1.00000 • 


2.00000 


12 


1.00000 


2.00000 


13 1 


1.00000 


2.00000 


14 


1.00000 : 


2.00000 


15 


1.00000 . 


2.00000 


16 


1.00000 


2.00000 


17 


1.00000 


2.00000 


18 


1.00000 


2.00000 


19 


1.00000 


2.00000 


20 


1.00000 


2.00000 



ejs 


* RT * 6 


7 


1 


. 99752 


. 99909 


2 


1 . 9801 7 


1.99179 


3 


2.91820 


2.96216 


4 


3 . 767 Q 0 


3.88039 


5 


4 . 48194 


4.70740 


6 


5.03626 


5 .40669 


7 


5.42996 


5 . 95698 


8 


5 . 6859 b 


6.35827 


9 


5.83874 


6.62918 


10 


5.92267 


6 .79868 


11 


5 . 96529 


6.89720 


12 


5.98538 


6.95055 


13 


5 . 99421 


6 .97755 


14 


5.99783 


6 . 99 o 36 


15 


5.99923 


6 . 996 ( j 6 


16 


5.99974 


6.99848 


17 


5 . 99992 


6 . 99944 


18 


5 . 99998 


6 . 99980 


19 


5 . 99999 


6 .99993 


20 


6. 00 0 00 


6 . 99998 



3 4 | 5 



. 95021 . 


. 98168 


.99326 


. 75106 


1 . 89011 


1.95263 


.32787 


2. 65200 


2.82818 


.68064 


3.21853 


3.56316 


.86538 


3 .58970 


4.12266 


.94930 


3.80457 


4.50670 


.98281 


3.91524 


4.74452 


.99471 


3.96637 


4.87789 


.99851 


' 3.98774 


4 .94598 


. 99962 


3 .99587 


4.97761 


.99991 


3.99871 


4.99151 


.99998 


3.99962 


4.99696 


•00000 i 


3 . 99990 


4 . 99898 


•00000 


3.99997 


4 . 99968 


•ooooo 


3.99999 


4.99990 


•00000 


4.00000 


4.99997 


•ooooo 


4.00000 


4 . 99999 


•ooooo 


4-00000 


5.00000 


•ooooo 


4.00000 


5.00000 


•ooooo 


4.00000 


5.00000 



8 


9 


10 


. 99966 


. 99988 


. 99995 


. 99665 


1 . 99864 


1 . 99946 


. 98269 


2 . 99241 


2 .99669 


. 94051 


3 . 9 7 i 1 6 


3.96635 


.84068 


4.91622 


4.95710 


.64964 


5.80053 


5.89001 


. 33627 


6.59375 


6 . 75987 


. 88331 


7 . 26965 


7.53965 


. 29076 


7.81420 


8.20683 


.57414 


8.22679 


8 . 74890 


. 75826 


8.52080 


9.16586 


.87017 


6.71779 


9.46908 


.93397 


8 . 84202 


9.67753 


.96815 


8 . 91587 


9 .81306 


. 96541 


8.95734 


9.89652 


.99364 


8 . 97937 


9.94526 


. 99736 


8 . 99048 


9.97230 


.99895 


8.99580 


9.98658 


. 99960 


8 . 99823 


9.99377 


.99986 


8.99928 


9 .99722 



1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

1 

2 

3 

4 

5 

6 

6 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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BJS 

1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

n 

12 

13 

14 

15 

19 

1 7 

ia 

19 

20 

BJS 

1 

2 

3 

4 

5 

a 

7 

8 

9 

10 

11 

12 

13 

14 

15 

19 

17 

18 

19 

20 



11 



RT 



12 13 14 



.99998 


. 99999 


1-00000 


1.00000 


1 


1.99976 

2 99657 

3 ! 99366 


1 .99991 


1 . 9999/ 


1 . 99999 


1 


2 . 9993v 


2 . 999/4 


2 . 99969 


2 


3.99710 


3 .99869 


3 . 99942 


3 


4*97655 


4.98950 


4 . 99495 


4 .99761 


4 


5.941Q3 


5 . 96916 


5.98422 


5.99206 


5 


6 .86242 


6 . 92334 


6 . 95633 


6.97765 


6 


7 .71923 


7 .83368 


7.90430 


7.94623 


i 7 


8 . 48724 


8 .67860 


8 .80455 


8 . 88418 


! 8 


9 . 1 4 6 7 3 


9 . 43641 


9 . 63873 


9.77476 


9 


9 . 68664 


10 .08916 


10.38/05 


10 - 59910 


10 


10 .10 75 7 


10 . 62759 


11.03369 


11.33906 


11 


10 .41888 


11 . 0 5 1 6 2 


11.57078 


11.96060 


12 


10.63 / 59 


11.37oo« 


11.997/4 


12 . 51615 


12 


10 • 78354 


11.59806 


12 . 32261 


12 .94572 


13 


10*87615 


11 . 75364 


12.55900 


13.27636 


13 


10-9320/ 


11 . 65493 


12.72351 


1 3 . 52044 


14 


10 .96426 


11.91/90 


12.83304 


13 . 69324 


14 


10 • 98195 


11 . 95532 


12.9028/ 


13.81059 


14 


10-99124 


11.97660 


12.94554 


13.68710 


14 



R r = i6 


17 


16 


19 




1.00000 


1-00000 


1-00000 


1.00000 


1 


2 . 0 0 0 0 U 


2.00000 


2.00000 


2.00000 


2 


2 . 99996 


2 . 9999V 


3.00000 


3.00000 


3 


3 . 999b9 


3 . 99995 


3 . 99998 


3 . 99999 


4 


4 . 99949 


4 . 99977 


4 . 99990 


4 .99995 


1 4 


5.99810 


5 . 99909 


5 . 9995 / 


5 .99960 


5 


6 . 99410 


6 . 997q3 


6 . 99853 


6 . 99928 


6 


7 . 98 4 1 u 


7 . 99160 


7 . 99563 


7.997// 


7 


8 . 96211 


8.97900 


8 . 98856 


8 . 99369 


8 


9.91661 


9.95286 


9 . 97320 


9.98504 


9, 


10.84142 


10.90375 : 


10 .94263 


10.96672 


10 . 


11 .71442 


11.81906 


11 . 88794 


11.9320 4 


11 , 


12.52130 


12 . 68406 


12. 79627 


12.87146 


12, 


13.24679 


13. 48319 


13 . 65368 


13. 77306 


13. 


13.87926 


14 . 20235 


14 .44560 


1 4 , 62333 


14. 


14.41252 


14.83090 


15.15895 


1 5 . 4 0 85 4 


15. 


14.84656 


15 . 36316 


15 . 78390 


16.11651 


16 . 


15.18721 


15.79914 


16.31525 


16.73814 


17. 


15 . 44486 


16 .14416 


16 . 75300 


17 . 26866 


17. 


15.63262 


16.40786 


17.10209 


1 7 . 7q8o5 


18. 



15 

0 00 00 

99999 
. 99999 
. 99974 
. 99889 
. 99610 
.98846 
.97046 
.93301 
.86316 
.74470 
. 55994 
.29233 
. 92912 
.46346 
. 89537 
. 23125 
. 48239 
. 66292 
.78770 
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